IDHwt glioblastomas can be stratified by their transcriptional response to standard treatment, with implications for targeted therapy.

BACKGROUND: Glioblastoma (GBM) brain tumors lacking IDH1 mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumors almost always fatally recur. RESULTS: Using RNA sequencing data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumors, we identify two responder subtypes based on longitudinal changes in gene expression. In two thirds of patients, a specific subset of genes is upregulated from primary to recurrence (Up responders), and in one third, the same genes are downregulated (Down responders), specifically in neoplastic cells. Characterization of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms that are associated with distinct changes in the tumor microenvironment. In Up responders, recurrent tumors are enriched in quiescent proneural GBM stem cells and differentiated neoplastic cells, with increased interaction with the surrounding normal brain and neurotransmitter signaling, whereas Down responders commonly undergo mesenchymal transition. ChIP-sequencing data from longitudinal GBM tumors suggests that the observed transcriptional reprogramming could be driven by Polycomb-based chromatin remodeling rather than DNA methylation. CONCLUSIONS: We show that the responder subtype is cancer-cell intrinsic, recapitulated in in vitro GBM cell models, and influenced by the presence of the tumor microenvironment. Stratifying GBM tumors by responder subtype may lead to more effective treatment.

The Epigenetic Evolution of Glioma Is Determined by the IDH1 Mutation Status and Treatment Regimen.

Tumor adaptation or selection is thought to underlie therapy resistance in glioma. To investigate longitudinal epigenetic evolution of gliomas in response to therapeutic pressure, we performed an epigenomic analysis of 132 matched initial and recurrent tumors from patients with IDH-wildtype (IDHwt) and IDH-mutant (IDHmut) glioma. IDHwt gliomas showed a stable epigenome over time with relatively low levels of global methylation. The epigenome of IDHmut gliomas showed initial high levels of genome-wide DNA methylation that was progressively reduced to levels similar to those of IDHwt tumors. Integration of epigenomics, gene expression, and functional genomics identified HOXD13 as a master regulator of IDHmut astrocytoma evolution. Furthermore, relapse of IDHmut tumors was accompanied by histologic progression that was associated with survival, as validated in an independent cohort. Finally, the initial cell composition of the tumor microenvironment varied between IDHwt and IDHmut tumors and changed differentially following treatment, suggesting increased neoangiogenesis and T-cell infiltration upon treatment of IDHmut gliomas. This study provides one of the largest cohorts of paired longitudinal glioma samples with epigenomic, transcriptomic, and genomic profiling and suggests that treatment of IDHmut glioma is associated with epigenomic evolution toward an IDHwt-like phenotype. SIGNIFICANCE: Standard treatments are related to loss of DNA methylation in IDHmut glioma, resulting in epigenetic activation of genes associated with tumor progression and alterations in the microenvironment that resemble treatment-naïve IDHwt glioma.

Driver Mutations Dictate the Immunologic Landscape and Response to Checkpoint Immunotherapy of Glioblastoma.

The composition of the tumor immune microenvironment (TIME) is considered a key determinant of patients' response to immunotherapy. The mechanisms underlying TIME formation and development over time are poorly understood. Glioblastoma (GBM) is a lethal primary brain cancer for which there are no curative treatments. GBMs are immunologically heterogeneous and impervious to checkpoint blockade immunotherapies. Utilizing clinically relevant genetic mouse models of GBM, we identified distinct immune landscapes associated with expression of EGFR wild-type and mutant EGFRvIII cancer driver mutations. Over time, accumulation of polymorphonuclear myeloid-derived suppressor cells (PMN-MDSC) was more pronounced in EGFRvIII-driven GBMs and was correlated with resistance to PD-1 and CTLA-4 combination checkpoint blockade immunotherapy. We determined that GBM-secreted CXCL1/2/3 and PMN-MDSC-expressed CXCR2 formed an axis regulating output of PMN-MDSCs from the bone marrow leading to systemic increase in these cells in the spleen and GBM tumor-draining lymph nodes. Pharmacologic targeting of this axis induced a systemic decrease in the numbers of PMN-MDSC, facilitated responses to PD-1 and CTLA-4 combination checkpoint blocking immunotherapy, and prolonged survival in mice bearing EGFRvIII-driven GBM. Our results uncover a relationship between cancer driver mutations, TIME composition, and sensitivity to checkpoint blockade in GBM and support the stratification of patients with GBM for checkpoint blockade therapy based on integrated genotypic and immunologic profiles.

GBMdeconvoluteR accurately infers proportions of neoplastic and immune cell populations from bulk glioblastoma transcriptomics data.

BACKGROUND: Characterizing and quantifying cell types within glioblastoma (GBM) tumors at scale will facilitate a better understanding of the association between the cellular landscape and tumor phenotypes or clinical correlates. We aimed to develop a tool that deconvolutes immune and neoplastic cells within the GBM tumor microenvironment from bulk RNA sequencing data. METHODS: We developed an IDH wild-type (IDHwt) GBM-specific single immune cell reference consisting of B cells, T-cells, NK-cells, microglia, tumor associated macrophages, monocytes, mast and DC cells. We used this alongside an existing neoplastic single cell-type reference for astrocyte-like, oligodendrocyte- and neuronal progenitor-like and mesenchymal GBM cancer cells to create both marker and gene signature matrix-based deconvolution tools. We applied single-cell resolution imaging mass cytometry (IMC) to ten IDHwt GBM samples, five paired primary and recurrent tumors, to determine which deconvolution approach performed best. RESULTS: Marker-based deconvolution using GBM-tissue specific markers was most accurate for both immune cells and cancer cells, so we packaged this approach as GBMdeconvoluteR. We applied GBMdeconvoluteR to bulk GBM RNAseq data from The Cancer Genome Atlas and recapitulated recent findings from multi-omics single cell studies with regards associations between mesenchymal GBM cancer cells and both lymphoid and myeloid cells. Furthermore, we expanded upon this to show that these associations are stronger in patients with worse prognosis. CONCLUSIONS: GBMdeconvoluteR accurately quantifies immune and neoplastic cell proportions in IDHwt GBM bulk RNA sequencing data and is accessible here: https://gbmdeconvoluter.leeds.ac.uk.

Glioma progression is shaped by genetic evolution and microenvironment interactions.

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

Central nervous system immune interactome is a function of cancer lineage, tumor microenvironment, and STAT3 expression.

BACKGROUNDImmune cell profiling of primary and metastatic CNS tumors has been focused on the tumor, not the tumor microenvironment (TME), or has been analyzed via biopsies.METHODSEn bloc resections of gliomas (n = 10) and lung metastases (n = 10) were analyzed via tissue segmentation and high-dimension Opal 7-color multiplex imaging. Single-cell RNA analyses were used to infer immune cell functionality.RESULTSWithin gliomas, T cells were localized in the infiltrating edge and perivascular space of tumors, while residing mostly in the stroma of metastatic tumors. CD163+ macrophages were evident throughout the TME of metastatic tumors, whereas in gliomas, CD68+, CD11c+CD68+, and CD11c+CD68+CD163+ cell subtypes were commonly observed. In lung metastases, T cells interacted with CD163+ macrophages as dyads and clusters at the brain-tumor interface and within the tumor itself and as clusters within the necrotic core. In contrast, gliomas typically lacked dyad and cluster interactions, except for T cell CD68+ cell dyads within the tumor. Analysis of transcriptomic data in glioblastomas revealed that innate immune cells expressed both proinflammatory and immunosuppressive gene signatures.CONCLUSIONOur results show that immunosuppressive macrophages are abundant within the TME and that the immune cell interactome between cancer lineages is distinct. Further, these data provide information for evaluating the role of different immune cell populations in brain tumor growth and therapeutic responses.FUNDINGThis study was supported by the NIH (NS120547), a Developmental research project award (P50CA221747), ReMission Alliance, institutional funding from Northwestern University and the Lurie Comprehensive Cancer Center, and gifts from the Mosky family and Perry McKay. Performed in the Flow Cytometry & Cellular Imaging Core Facility at MD Anderson Cancer Center, this study received support in part from the NIH (CA016672) and the National Cancer Institute (NCI) Research Specialist award 1 (R50 CA243707). Additional support was provided by CCSG Bioinformatics Shared Resource 5 (P30 CA046592), a gift from Agilent Technologies, a Research Scholar Grant from the American Cancer Society (RSG-16-005-01), a Precision Health Investigator Award from University of Michigan (U-M) Precision Health, the NCI (R37-CA214955), startup institutional research funds from U-M, and a Biomedical Informatics & Data Science Training Grant (T32GM141746).

Single-cell multimodal glioma analyses identify epigenetic regulators of cellular plasticity and environmental stress response.

Glioma intratumoral heterogeneity enables adaptation to challenging microenvironments and contributes to therapeutic resistance. We integrated 914 single-cell DNA methylomes, 55,284 single-cell transcriptomes and bulk multi-omic profiles across 11 adult IDH mutant or IDH wild-type gliomas to delineate sources of intratumoral heterogeneity. We showed that local DNA methylation disorder is associated with cell-cell DNA methylation differences, is elevated in more aggressive tumors, links with transcriptional disruption and is altered during the environmental stress response. Glioma cells under in vitro hypoxic and irradiation stress increased local DNA methylation disorder and shifted cell states. We identified a positive association between genetic and epigenetic instability that was supported in bulk longitudinally collected DNA methylation data. Increased DNA methylation disorder associated with accelerated disease progression and recurrently selected DNA methylation changes were enriched for environmental stress response pathways. Our work identified an epigenetically facilitated adaptive stress response process and highlights the importance of epigenetic heterogeneity in shaping therapeutic outcomes.

Radiotherapy is associated with a deletion signature that contributes to poor outcomes in patients with cancer.

Ionizing radiation causes DNA damage and is a mainstay for cancer treatment, but understanding of its genomic impact is limited. We analyzed mutational spectra following radiotherapy in 190 paired primary and recurrent gliomas from the Glioma Longitudinal Analysis Consortium and 3,693 post-treatment metastatic tumors from the Hartwig Medical Foundation. We identified radiotherapy-associated significant increases in the burden of small deletions (5-15 bp) and large deletions (20+ bp to chromosome-arm length). Small deletions were characterized by a larger span size, lacking breakpoint microhomology and were genomically more dispersed when compared to pre-existing deletions and deletions in non-irradiated tumors. Mutational signature analysis implicated classical non-homologous end-joining-mediated DNA damage repair and APOBEC mutagenesis following radiotherapy. A high radiation-associated deletion burden was associated with worse clinical outcomes, suggesting that effective repair of radiation-induced DNA damage is detrimental to patient survival. These results may be leveraged to predict sensitivity to radiation therapy in recurrent cancer.

Perspective of mesenchymal transformation in glioblastoma.

Despite aggressive multimodal treatment, glioblastoma (GBM), a grade IV primary brain tumor, still portends a poor prognosis with a median overall survival of 12-16 months. The complexity of GBM treatment mainly lies in the inter- and intra-tumoral heterogeneity, which largely contributes to the treatment-refractory and recurrent nature of GBM. By paving the road towards the development of personalized medicine for GBM patients, the cancer genome atlas classification scheme of GBM into distinct transcriptional subtypes has been considered an invaluable approach to overcoming this heterogeneity. Among the identified transcriptional subtypes, the mesenchymal subtype has been found associated with more aggressive, invasive, angiogenic, hypoxic, necrotic, inflammatory, and multitherapy-resistant features than other transcriptional subtypes. Accordingly, mesenchymal GBM patients were found to exhibit worse prognosis than other subtypes when patients with high transcriptional heterogeneity were excluded. Furthermore, identification of the master mesenchymal regulators and their downstream signaling pathways has not only increased our understanding of the complex regulatory transcriptional networks of mesenchymal GBM, but also has generated a list of potent inhibitors for clinical trials. Importantly, the mesenchymal transition of GBM has been found to be tightly associated with treatment-induced phenotypic changes in recurrence. Together, these findings indicate that elucidating the governing and plastic transcriptomic natures of mesenchymal GBM is critical in order to develop novel and selective therapeutic strategies that can improve both patient care and clinical outcomes. Thus, the focus of our review will be on the recent advances in the understanding of the transcriptome of mesenchymal GBM and discuss microenvironmental, metabolic, and treatment-related factors as critical components through which the mesenchymal signature may be acquired. We also take into consideration the transcriptomic plasticity of GBM to discuss the future perspectives in employing selective therapeutic strategies against mesenchymal GBM.

Very low mutation burden is a feature of inflamed recurrent glioblastomas responsive to cancer immunotherapy.

Several immunotherapy clinical trials in recurrent glioblastoma have reported long-term survival benefits in 10-20% of patients. Here we perform genomic analysis of tumor tissue from recurrent WHO grade IV glioblastoma patients acquired prior to immunotherapy intervention. We report that very low tumor mutation burden is associated with longer survival after recombinant polio virotherapy or after immune checkpoint blockade in recurrent glioblastoma patients. A relationship between tumor mutation burden and survival is not observed in cohorts of immunotherapy naïve newly diagnosed or recurrent glioblastoma patients. Transcriptomic analyses reveal an inverse relationship between tumor mutation burden and enrichment of inflammatory gene signatures in cohorts of recurrent, but not newly diagnosed glioblastoma tumors, implying that a relationship between tumor mutation burden and tumor-intrinsic inflammation evolves upon recurrence.

Whole transcriptome signature for prognostic prediction (WTSPP): application of whole transcriptome signature for prognostic prediction in cancer.

Developing prognostic biomarkers for specific cancer types that accurately predict patient survival is increasingly important in clinical research and practice. Despite the enormous potential of prognostic signatures, proposed models have found limited implementations in routine clinical practice. Herein, we propose a generic, RNA sequencing platform independent, statistical framework named whole transcriptome signature for prognostic prediction to generate prognostic gene signatures. Using ovarian cancer and lung adenocarcinoma as examples, we provide evidence that our prognostic signatures overperform previous reported signatures, capture prognostic features not explained by clinical variables, and expose biologically relevant prognostic pathways, including those involved in the immune system and cell cycle. Our approach demonstrates a robust method for developing prognostic gene expression signatures. In conclusion, our statistical framework can be generally applied to all cancer types for prognostic prediction and might be extended to other human diseases. The proposed method is implemented as an R package (PanCancerSig) and is freely available on GitHub ( https://github.com/Cheng-Lab-GitHub/PanCancer_Signature ).

Comparative Molecular Life History of Spontaneous Canine and Human Gliomas.

Sporadic gliomas in companion dogs provide a window on the interaction between tumorigenic mechanisms and host environment. We compared the molecular profiles of canine gliomas with those of human pediatric and adult gliomas to characterize evolutionarily conserved mammalian mutational processes in gliomagenesis. Employing whole-genome, exome, transcriptome, and methylation sequencing of 83 canine gliomas, we found alterations shared between canine and human gliomas such as the receptor tyrosine kinases, TP53 and cell-cycle pathways, and IDH1 R132. Canine gliomas showed high similarity with human pediatric gliomas per robust aneuploidy, mutational rates, relative timing of mutations, and DNA-methylation patterns. Our cross-species comparative genomic analysis provides unique insights into glioma etiology and the chronology of glioma-causing somatic alterations.

Multifactorial Deep Learning Reveals Pan-Cancer Genomic Tumor Clusters with Distinct Immunogenomic Landscape and Response to Immunotherapy.

PURPOSE: Tumor genomic features have been of particular interest because of their potential impact on the tumor immune microenvironment and response to immunotherapy. Due to the substantial heterogeneity, an integrative approach incorporating diverse molecular features is needed to characterize immunologic features underlying primary resistance to immunotherapy and for the establishment of novel predictive biomarkers. EXPERIMENTAL DESIGN: We developed a pan-cancer deep machine learning model integrating tumor mutation burden, microsatellite instability, and somatic copy-number alterations to classify tumors of different types into different genomic clusters, and assessed the immune microenvironment in each genomic cluster and the association of each genomic cluster with response to immunotherapy. RESULTS: Our model grouped 8,646 tumors of 29 cancer types from The Cancer Genome Atlas into four genomic clusters. Analysis of RNA-sequencing data revealed distinct immune microenvironment in tumors of each genomic class. Furthermore, applying this model to tumors from two melanoma immunotherapy clinical cohorts demonstrated that patients with melanoma of different genomic classes achieved different benefit from immunotherapy. Interestingly, tumors in cluster 4 demonstrated a cold immune microenvironment and lack of benefit from immunotherapy despite high microsatellite instability burden. CONCLUSIONS: Our study provides a proof for principle that deep learning modeling may have the potential to discover intrinsic statistical cross-modality correlations of multifactorial input data to dissect the molecular mechanisms underlying primary resistance to immunotherapy, which likely involves multiple factors from both the tumor and host at different molecular levels.

Longitudinal molecular trajectories of diffuse glioma in adults.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

Single-cell RNA sequencing reveals the impact of chromosomal instability on glioblastoma cancer stem cells.

BACKGROUND: Intra-tumor heterogeneity stems from genetic, epigenetic, functional, and environmental differences among tumor cells. A major source of genetic heterogeneity comes from DNA sequence differences and/or whole chromosome and focal copy number variations (CNVs). Whole chromosome CNVs are caused by chromosomal instability (CIN) that is defined by a persistently high rate of chromosome mis-segregation. Accordingly, CIN causes constantly changing karyotypes that result in extensive cell-to-cell genetic heterogeneity. How the genetic heterogeneity caused by CIN influences gene expression in individual cells remains unknown. METHODS: We performed single-cell RNA sequencing on a chromosomally unstable glioblastoma cancer stem cell (CSC) line and a control normal, diploid neural stem cell (NSC) line to investigate the impact of CNV due to CIN on gene expression. From the gene expression data, we computationally inferred large-scale CNVs in single cells. Also, we performed copy number adjusted differential gene expression analysis between NSCs and glioblastoma CSCs to identify copy number dependent and independent differentially expressed genes. RESULTS: Here, we demonstrate that gene expression across large genomic regions scales proportionally to whole chromosome copy number in chromosomally unstable CSCs. Also, we show that the differential expression of most genes between normal NSCs and glioblastoma CSCs is largely accounted for by copy number alterations. However, we identify 269 genes whose differential expression in glioblastoma CSCs relative to normal NSCs is independent of copy number. Moreover, a gene signature derived from the subset of genes that are differential expressed independent of copy number in glioblastoma CSCs correlates with tumor grade and is prognostic for patient survival. CONCLUSIONS: These results demonstrate that CIN is directly responsible for gene expression changes and contributes to both genetic and transcriptional heterogeneity among glioblastoma CSCs. These results also demonstrate that the expression of some genes is buffered against changes in copy number, thus preserving some consistency in gene expression levels from cell-to-cell despite the continuous change in karyotype driven by CIN. Importantly, a gene signature derived from the subset of genes whose expression is buffered against copy number alterations correlates with tumor grade and is prognostic for patient survival that could facilitate patient diagnosis and treatment.

Hypoxia-Induced VISTA Promotes the Suppressive Function of Myeloid-Derived Suppressor Cells in the Tumor Microenvironment.

Tumor hypoxia is a negative prognostic factor that is implicated in oncogenic signal activation, immune escape, and resistance to treatment. Identifying the mechanistic role of hypoxia in immune escape and resistance to immune-checkpoint inhibitors may aid the identification of therapeutic targets. We and others have shown that V-domain Ig suppressor of T-cell activation (VISTA), a negative checkpoint regulator in the B7 family, is highly expressed in the tumor microenvironment in tumor models and primary human cancers. In this study, we show that VISTA and HIF1α activity are correlated in a cohort of colorectal cancer patients. High VISTA expression was associated with worse overall survival. We used the CT26 colon cancer model to investigate the regulation of VISTA by hypoxia. Compared with less hypoxic tumor regions or draining lymph nodes, regions of profound hypoxia in the tumor microenvironment were associated with increased VISTA expression on tumor-infiltrating myeloid-derived suppressor cells (MDSC). Using chromatin immunoprecipitation and genetic silencing, we show that hypoxia-inducible factor (HIF)-1α binding to a conserved hypoxia response element in the VISTA promoter upregulated VISTA on myeloid cells. Further, antibody targeting or genetic ablation of VISTA under hypoxia relieved MDSC-mediated T-cell suppression, revealing VISTA as a mediator of MDSC function. Collectively, these data suggest that targeting VISTA may mitigate the deleterious effects of hypoxia on antitumor immunity.

Cytokine sensitivity screening highlights BMP4 pathway signaling as a therapeutic opportunity in ER+ breast cancer.

Despite the success of approved systemic therapies for estrogen receptor α (ER)-positive breast cancer, drug resistance remains common. We hypothesized that secreted factors from the human tumor microenvironment could modulate drug resistance. We previously screened a library of 297 recombinant-secreted microenvironmental proteins for the ability to confer resistance to the anti-estrogen fulvestrant in 2 ER+ breast cancer cell lines. Herein, we considered whether factors that enhanced drug sensitivity could be repurposed as therapeutics and provide leads for drug development. Screening data revealed bone morphogenic protein (BMP)4 as a factor that inhibited cell growth and synergized with approved anti-estrogens and cyclin-dependent kinase 4/6 inhibitors (CDK4/6i). BMP4-mediated growth inhibition was dependent on type I receptor activin receptor-like kinase (ALK)3-dependent phosphorylation (P) of mothers against decapentaplegic homolog (SMAD/P-SMAD)1 and 5, which could be reversed by BMP receptor inhibitors and ALK3 knockdown. The primary effect of BMP4 on cell fate was cell-cycle arrest, in which RNA sequencing, immunoblot analysis, and RNA interference revealed to be dependent on p21WAF1/Cip1 upregulation. BMP4 also enhanced sensitivity to approved inhibitors of mammalian target of rapamycin complex 1 and CDK4/6 via ALK3-mediated P-SMAD1/5 and p21 upregulation in anti-estrogen-resistant cells. Patients bearing primary ER+ breast tumors, exhibiting a transcriptomic signature of BMP4 signaling, had improved disease outcome following adjuvant treatment with anti-estrogen therapy, independently of age, tumor grade, and tumor stage. Furthermore, a transcriptomic signature of BMP4 signaling was predictive of an improved biologic response to the CDK4/6i palbociclib, in combination with an aromatase inhibitor in primary tumors. These findings highlight BMP4 and its downstream pathway activation as a therapeutic opportunity in ER+ breast cancer.-Shee, K., Jiang, A., Varn, F. S., Liu, S., Traphagen, N. A., Owens, P., Ma, C. X., Hoog, J., Cheng, C., Golub, T. R., Straussman, R., Miller, T. W. Cytokine sensitivity screening highlights BMP4 pathway signaling as a therapeutic opportunity in ER+ breast cancer.

A B cell-derived gene expression signature associates with an immunologically active tumor microenvironment and response to immune checkpoint blockade therapy.

Immune checkpoint inhibitors have shown great potential in treating solid tumors, inducing durable remission and prolonged survival time in responders. Despite their promise, a large fraction of patients remains unresponsive to these treatments highlighting the need for biomarkers that can predict patient sensitivity. Pre-treatment gene expression profiles for patients receiving immune checkpoint inhibitors have recently become available, establishing a new medium by which to discover biomarkers that predict therapy response. In this study, we mined for transcriptomic correlates of response by applying immune cell-derived gene expression signatures to publicly available datasets containing matched gene expression and response efficacy information. These datasets were comprised of urothelial carcinoma patients receiving anti-PD-L1 (n = 25), melanoma patients receiving anti-PD-1 (n = 28), and melanoma patients receiving anti-CTLA-4 (n = 42). We identified one signature, derived from a subpopulation of B cells, with scores that were significantly and reproducibly elevated in patients experiencing clinical benefit following therapy targeting the PD-1/PD-L1 axis and were additionally elevated in patients responsive to anti-CTLA-4 therapy. Multivariate models revealed that this signature was associated with response independent of other response-predictive biomarkers, including tumor mutation burden. Functional annotation of the signature revealed it to be associated with features indicative of an immunologically active microenvironment, including B and T cell activation as well as antigen presentation activity. The preliminary findings presented detail a transcriptomic signature associated with response to multiple checkpoint inhibitors and suggest novel biological associations that warrant further investigation.

Genomic Characterization of Six Virus-Associated Cancers Identifies Changes in the Tumor Immune Microenvironment and Altered Genetic Programs.

Viruses affect approximately 20% of all human cancers and induce expression of immunogenic viral oncoproteins that make these tumors potent targets for immune checkpoint inhibitors. In this study, we apply computational tools to The Cancer Genome Atlas (TCGA) and other genomic datasets to define how virus infection shapes the tumor immune microenvironment and genetic architecture of 6 virus-associated tumor types. Across cancers, the cellular composition of the microenvironment varied by viral status, with virus-positive tumors often exhibiting increased infiltration of cytolytic cell types compared with their virus-negative counterparts. Analyses of the infiltrating T-cell receptor repertoire in these patients revealed that Epstein-Barr virus infection was associated with decreased receptor diversity in multiple cancers, suggesting an antigen-driven clonal T-cell response. Tissue-specific gene-expression signatures capturing virus-associated transcriptomic changes successfully predicted virus status in independent datasets and were associated with both immune- and proliferation-related features that were predictive of patient prognosis. Together, the analyses presented suggest viruses have distinct effects in different tumors, with implications for immunotherapy.Significance: This study utilizes TCGA and other genomic datasets to further our understanding of how viruses affect the tumor immune response in different cancer types.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/78/22/6413/F1.large.jpg Cancer Res; 78(22); 6413-23. ©2018 AACR.

Computational immune profiling in lung adenocarcinoma reveals reproducible prognostic associations with implications for immunotherapy.

Non-small cell lung cancer is one of the leading causes of cancer-related death in the world. Lung adenocarcinoma, the most common type of non-small cell lung cancer, has been well characterized as having a dense lymphocytic infiltrate, suggesting that the immune system plays an active role in shaping this cancer's growth and development. Despite these findings, our understanding of how this infiltrate affects patient prognosis and its association with lung adenocarcinoma-specific clinical factors remains limited. To address these questions, we inferred the infiltration level of six distinct immune cell types from a series of four lung adenocarcinoma gene expression datasets. We found that naive B cell, CD8+ T cell, and myeloid cell-derived expression signals of immune infiltration were significantly predictive of patient survival in multiple independent datasets, with B cell and CD8+ T cell infiltration associated with prolonged prognosis and myeloid cell infiltration associated with shorter survival. These associations remained significant even after accounting for additional clinical variables. Patients stratified by smoking status exhibited decreased CD8+ T cell infiltration and altered prognostic associations, suggesting potential immunosuppressive mechanisms in smokers. Survival analyses accounting for immune checkpoint gene expression and cellular immune infiltrate indicated checkpoint protein-specific modulatory effects on CD8+ T cell and B cell function that may be associated with patient sensitivity to immunotherapy. Together, these analyses identified reproducible associations that can be used to better characterize the role of immune infiltration in lung adenocarcinoma and demonstrate the utility in using computational approaches to systematically characterize tissue-specific tumor-immune interactions.